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Abstract. Ingrid Carbone introduced the notion of so-called LS-sequences of 
points, which are obtained by a generalization of Kakutani's interval splitting 
procedure. Under an appropriate choice of the parameters L and S, such se- 
quences have low discrepancy, which means that they are natural candidates 
for Quasi-Montc Carlo integration. It is tempting to assume that LS-sequences 
can be combined coordinatewise to obtain a multidimensional low-discrepancy 
sequence. However, in the present paper we prove that this is not always 
the case: if the parameters L±,Si and L2,S I 2 of two one-dimensional low- 
discrepancy LS-sequences satisfy certain number-theoretic conditions, then 
their two-dimensional combination is not even dense in [0, l] 2 . 



1. Introduction and statement of results 

For two points a, b <E [0, l) d we write a < b and a < b if the correspond- 
ing inequalities hold in each coordinate; furthermore, we write [a, b) for the set 
{x £ [0,l) d : a < x < b}, and call such a set a (d-dimensional) interval. 
We denote by 1/ the indicator function of the set / C [0, l) d and by the d- 
dimensional Lebesgue measure. We will sometimes write for the ci-dimensional 
vector (0, . . . , 0). 

A sequence (i n ) n gn of points in [0,l] d is called uniformly distributed modulo 1 
(u.d. mod 1) if 

J™ AT = A d(M)) 

for all d-dimensional intervals [a, b) C [0, l) d . A further characterization of u.d. is 
given by the following well-known result of Weyl |13j : a sequence (x n ) n ^jq of points 
in [0, l) d is u.d. mod 1 if and only if for every continuous function / on [0, l) d the 
relation 

li m En=l -f( Xn ^ = f f( x )dx 

holds. Although this theorem shows the possibility of using u.d. point sequences 
for numerical integration — a method usually called Quasi-Monte Carlo (QMC) 
integration — it does not give any information on the integration error. 
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The Koksma-Hlawka inequality [8] states that this integration error can be 
bounded by the product of the variation of / (in the sense of Hardy and Krause), 
denoted by V(f), and the so-called star-discrepancy D* N of the point sequence 
(x„)„ eN . More precisely, 



The Koksma-Hlawka inequality suggest that for QMC integration one should use 
a sequence of points whose discrepancy is as small as possible. The best known 
point sequences achieve a discrepancy of order 0(N^ 1 (log N) d ); such sequences are 
called low- discrepancy sequences. Note that this convergence rate is for all d > 1 
better than that of the probabilistic error of Monte Carlo integration, where a se- 
quence of random points is used instead of a deterministic one. QMC integration 
can be successfully applied in several different areas of applied mathematics, for 
example in actuarial or financial mathematics, where frequently high-dimensional 
integration problems arise (see e.g. [21 HI])- For more information on discrepancy 
theory, low-discrepancy sequences and QMC integration see [3 [10] . 

A popular approach to construct <i-dimensional low-discrepancy sequences is 
to combine d one-dimensional low-discrepancy sequences; this works for example 
for the so-called Halton sequence, which is obtained by joining one-dimensional 
van der Corput sequences coordinatewise. In the present paper we show that 
this construction principle is not generally applicable for a special class of one- 
dimensional low-discrepancy sequence, so-called LS-sequences. We prove that the 
limit distribution of a multidimensional LS-sequences (composed coordinatewise 
from one-dimensional low-discrepancy LS-sequences) can spectacularly fail to be 
u.d., if there is a certain number-theoretic connection between the parameters of 
the one-dimensional sequences. To explain the construction of LS-sequences we 
need some definitions. 

Definition 1 (Kakutani splitting procedure). If a £ (0, 1) and ir — {[£j_x,ij) : 1 < 
i < k} is any partition of [0, 1), then om denotes its so-called a-refinement, which 
is obtained by subdividing all intervals of ir having maximal length into two parts, 
proportional to a and 1 — a, respectively. The so-called Kakutani 's sequence of par- 
titions (a™w)„ e N is obtained as the successive a-refinement of the trivial partition 



The notion of a-refinements can be generalized in a natural way to so-called 
p-refinements. 

Definition 2 (p-refinement). Let p denote a non-trivial finite partition of [0, 1). 
Then the p-refinement of a partition tt of [0, 1), denoted by pir, is given by subdi- 
viding all intervals of maximal length positively homothetically to p. Note that the 
a-refinement is a special case with p = {[0, a), [a, 1)}. 




where is given by 



D% = Dx(xi,...,x N ) = sup 

a£[0,l) d 




N 



c = {[0,l)}. 
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By a classical result of Kakutani [9] , for any a the sequence of partitions (a"w) n eN 
is uniformly distributed, which means that for every interval [a, b] C [0, 1], 

lim — — — = o — a, 

n->oo fe(n) 

where fc(n) denotes the number of intervals in a"u = {[i"_ l7 i"), 1 < i < fc(n)}. 
The same result holds for any sequence of p- refinements of u>, due to a result of 
Volcic [12] (see also [HE]). A multidimensional generalization of p-refinements has 
been introduced by Carbone and Volcic [5]. A special case of a p- refinement is the 
so-called LS-sequence of partitions. This sequence of partitions has been introduced 
by Carbone [3]. 

Definition 3 (LS-sequence of partitions). An LS-sequence of partitions 
(p£ gu)neN is the successive p-refinement of the trivial partition uj , where pl,s con- 
sists of L + S intervals such that the first L > intervals of pl,s have length (3 and 
the remaining S > intervals have length /3 2 . Note that necessarily L/3 + S/3 2 = 1 
holds, and consequently for each pair [L, S) of parameters there exists exactly one 
corresponding number (3. 

It can easily be seen that for every n the partition p\ s uj consists only of intervals 
having either length j3 n or f3 n+1 . This fact makes the analysis of LS-sequences 
relatively simple, in comparison to the analysis of general p-refinements. We denote 
by t n the total number of intervals of p£ s w ' anc ^ correspondingly let l n and s n be 
the number of long and short intervals after n steps, respectively (more precisely, 
l n is the number of intervals of length f3 n , and s n is the number of intervals of 
length f3 n+1 ). It is easy to see that these sequences satisfy the following recurrence 
relations (see [4]): 



t n 


= Lt n -i 


+ St n -2, 


I n 


= Ll n -\ 


■f Sl n -2, 


Sn 


= Ls n -i 


+ 


for n > 2, where t\ = L + S, to 


= Mi = 


L, Iq, = 



these binary recurrences we obtain explicit formulas for t n , l n and s T . 

(1) tn = T f3- n + T 1 (Sf3) n , T = 



L + 2S+ VU+4S 



2y/L 2 + AS 



n 

(2) l n = X /3- n + X 1 {-Sf3) n , A 

Ai 

(3) s n = a f3- n + a 1 (-Sf3) n , a Q 

o-i 



-L-2S+\/L 2 +4S 
2y/L 2 + 45 ' 
L + \/L 2 + 45 
2VL 2 +45~' 



-L + ^L 2 + 4S 



2VL 2 + iS 



2S+VL 2 +4S 
2VL 2 +45 1- ' 



-2S + VL 2 +45' 
2^L 2 + AS 
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We can generate a sequence of points from a sequence of partitions by ordering 
the left endpoints of the intervals in the partition. The following rule by Carbone [4] 
defines the so-called LS-sequence of points. 

Definition 4 (LS-sequence of points). Given an LS-sequence of partitions 
(p2 s<^)neN we define the corresponding LS-sequence of points (££ s )neN as follows: 
the first t\ points are the left endpoints of the partition pl,s^ ordered by magnitude. 
We denote this ordered set of points by s . 

For n > 1 we define A^'g 1 = 5, • ■ • , CiTj 1 } inductively as the ordered set of 
the left endpoints of the intervals of p£ 5^ in the following way: 

A»+i _/<r0 tt„-l 

v^ 1 (tis) , - • • , (d^ 1 ) , . . . , ^to 1 , • • • , ^o 1 (d?/) , 

• • • • . & Ws 1 ) , • • • , « , . ■ • , tfft 1 (W) >> 

w/iere 

^(a?) =x + i/3"+j/3 n+1 . 

For more details on the definition of LS-sequences of points, and on the properties 
of such sequences, see [3j S] . 

Next, we recall the definition of the well-known van der Corput sequence in base 
b > 2, b £ N. For every n € No, the unique digit expansion of n in base b is given 

by 

i>0 

where n% € {0, 1, . . . , b— 1}, i > 0. For ngNo we define the radical-inverse function 
(or Monna map) 4>b{n) : No — > [0, 1) by 

(4) 06 (n) = 06 I XI ] := S nift - * -1 . 

\i>0 / i>0 

We call x a 6-adic rational if x = ab~ c , where a and c are positive integers and 
< a < b c . Note that <f>b(n) maps N onto the &-adic rationals in [0, 1), and therefore 
the image of N under <fib(n) is dense in [0, 1). 

Definition 5. The van der Corput sequence in base b is defined as (4>b(ri))neN- 

Note that the definition of the van der Corput sequence in base b > 2 coincides 
with the definition of the LS-sequence of points with parameters L — b and 5 = 0. 
Thus LS-sequences can be seen as a generalization of the van der Corput sequence. 
A remarkable property of van der Corput sequences is, that several van der Corput 
sequences in pairwise coprime bases can be combined coordinatewise to a multi- 
dimensional sequence, the so-called Halton sequence, which is a low-discrepancy 
sequence. As mentioned above this means that the discrepancy of a Halton se- 
quence is of asymptotic order 0(N~ 1 (log N) d ), where N is the number of points 
and d denotes the dimension, which together with the Koksma-Hlawka inequality 
makes it a perfect candidate for Quasi-Monte Carlo integration (for details on the 
properties of van der Corput and Halton sequences, see [71 HO]). 
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If we assume S > 1, then by a result of Carbone [4,, one-dimensional LS- 
sequences are low-discrepancy sequences if and only if L > 5—1. Thus it is 
tempting to assume that several LS-sequences can be combined coordinatewise in 
order to obtain multidimensional low-discrepancy sequences. If this was be the 
case, then this method would produce a new parametric class of multidimensional 
low-discrepancy sequences. However, even in the case of the combination of van 
der Corput sequences (which are a special case of LS-sequences, as mentioned be- 
fore), the bases b\, . . . , bd cannot be chosen arbitrarily, but have to satisfy a certain 
number-theoretic condition (they have to be coprime). A similar restriction can be 
expected in the case of combining LS-sequences. 

In a talk in Graz in June 2012, Maria Rita Iaco presented several numerical 
examples of the asymptotic distribution of two-dimensional LS-sequences. In some 
cases they showed "random" behavior, while in others (for example when combining 
the sequence with parameters (1,1) and the sequence with parameters (4,1)) the dis- 
tribution seemed to be rather erratic. Obviously the reason for this behavior is that 
there is a multiplicative relation between the solutions of the equations x + x 2 = 1 
and Ax + x 2 = 1, which define the lengths of the intervals for the LS-sequences with 
parameters (1,1) and (4,1), respectively. The purpose of the present paper is to 
prove that in fact the two-dimensional LS-sequences is not uniformly distributed 
(and not even dense) in [0, l] 2 if such a multiplicative relation exists. Further- 
more, in a second theorem we show that if the parameters of two one-dimensional 
LS-sequences have a greatest common divisor (gcd) which is greater than 1, then 
the resulting two-dimensional LS-sequence is also not dense in [0, l] 2 . This second 
result generalizes the requirement of having coprime bases of the van der Corput 
sequences, in order to obtain a low-discrepancy Halton sequence by joining them 
coordinatewise . 

The formal definition of a multidimensional LS-sequence can be given as follows. 

Definition 6. Let B — ((Li, Si), . . . , (Lj, S^f) be an ordered d-tuple of pairs, 
(Li, Si) such that L t > 0, Si > and Li + Si > 2 for all i. Then we define 
the d- dimensional LS-sequence in base B as the sequence 

The following theorem states that a two-dimensional LS-sequences in bases B = 
((L\, Si), (L2, S2)), where the one-dimensional components are low-discrepancy se- 

quences, is not dense in [0, l] 2 if there exist integers m and k such that ^L+i € Q- 

For example, in the case {L\, Si) — (1, 1) and (1/2,52) = (4, 1) we have /?2 = P\- 

Theorem 1. Let B — ((£1, Si), (L2, S2)) with Li > Si — 1 > and assume that 

there exist integers m and k such that * + i G Q. Then the two-dimensional LS- 

"2 

sequence £g is not uniformly distributed, and not even dense in [0, l] 2 . 

On the other hand, we have not been able to derive any positive results, proving 
uniform distribution of a LS-sequence for an appropriate choice of Li, S\ and L2, S2 
(except for the case of the Halton sequence). So up to date not a single example 
of parameters Li,Si,L 2 , S2 is known, for which either Si ^ or S2 7^ and the 
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corresponding two-dimensional LS-sequence is uniformly distributed. 

Note that Theorem [T] can also be applied to the multidimensional case, since for 
any multidimensional sequence of points, which is uniformly distributed, all lower- 
dimensional projections also have to be uniformly distributed. More precisely, we 
immediately get the following corollary. 

Corollary 1. Let B = ((£1, Si), . . . , (£<£, Sd)) with L ; > S^ — 1 > and assume that 

flfe+i 

there exist numbers u,w £ {1, . . . , d} and integers m and k such that ^m+i G Q- 
Then the d-dimensional LS-sequence Q is not uniform distributed, and not even 
dense in [0, l] d . 

The next theorem characterizes another class of two-dimensional LS-sequences, 
which are not dense in [0, l] 2 . 

Theorem 2. Let B — ((L\, Si), (L2, S2)) and assume that gcd(Li, Si, L2, S2) > 1- 
Then the two-dimensional LS-sequence £g is not dense in [0, l] 2 . 

Note that Theorem [5] also includes the case of the Halton sequence. As above 
we can state a corollary which describes the d-dimensional situation. 

Corollary 2. Let B = ((£1, Si), . . . , (Ld, Sd)) and assume that gcd(L\, Si, Lj, Sj) > 
1 for some i,j 6 {l,...,d}, i ^ j '. Then the d-dimensional LS-sequence £g is not 
dense in [0, l] d . 

In the next section we provide several auxiliary results concerning one-dimensional 
LS-sequences. These lemmas will be essential in the proofs of Theorem [1] and The- 
orem [21 which are presented in Section [3] 

2. Points in elementary intervals 

Before we define elementary intervals and prove some of their properties, we 
prove the following recurrence relations for the sequences t n and l n . 

Lemma 1. We have 

t n +i = t n + (L + S - l)l n and l n+ i =*„ + (£- l)l„ 

for all n > 0. 

Proof. We prove this claim by induction on n. For n = we have 
t 1 =L + S=l + (L + S-l) = t + (L + S- l)l 

and 

h = L = 1 + (L - 1) = t + (L - l)l . 
Assume the claim is true for all indices < n. Then we have 

t n+1 = Ltn + St n -i = L(t n -i + {L + S- l)l n -l) + S(t„_ 2 + {L + S- l)l n -2) = 
Ltn-i + S<„_ 2 + (L + S- l){Ll n - X + Sl n -2) =t n + (L + S- l)l n 

and 

ln+i = Ll n + Sl n -x = L(t n -i + (L — l)Z n _i) + S(t n —2 + (L — l)l n -2) = 

Lt^ + Si„_ 2 + (L— l)(Ll n -i + Sl„- 2 ) = t n + {L- l)l n 

which proves the lemma. □ 
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We will also need a lemma on the irrationality of ft. 

Lemma 2. Let L > 5 — 1 > 0, then f3 k is irrational for every positive integer k. 

Proof. First, we prove that (3 is irrational. In particular we have to prove that 
L 2 + 45 is not a square. Therefore we note that 

I? < L 2 + 45 < L 2 + AL + 4 = (L + 2) 2 , 

thus L 2 + 45 = (L + l) 2 , i.e. 5 = 2 -^±I g Z. Hence f3 is irrational. 

Now suppose that f3 k is a rational for some k > 1. Let us assume that k is 
minimal. Then j3 is an algebraic integer of degree fc, which is a contradiction unless 
k = 2, since /? is the solution of a quadratic equation. But k = 2 yields 

which is not possible unless L — 0. □ 
We call an interval elementary, if it is an element of p£ s u; for some n. Equiv- 

(k) 

alently we can define elementary intervals as all intervals of the form 
Kl.S'£l.s + P k ) f° r some k, where x < Ik- If [££ s>£l.s + fi k ) 1S an elementary 
interval, then there necessarily exists an integer y < t^ such that s + /3 = t; y L s . 
Note that in the case of the van der Corput sequence the elementary intervals de- 
fined in this way coincide with the usual ones. 

In order to count points in elementary intervals we need a method to decide 
whether a point ^ s is contained in some given elementary interval or not. In 
the case of the van der Corput sequence this can be achieved by considering digit 
expansions as in ^j, which motivates the following construction. Let N > be a 
fixed integer and let n be such that t n < N < t n+ i. We construct two sequences 
(efc)o<fc<n and (r]k)o<k<n recursively in the following way: We put 

N n - U 



N n — N, e n = 1, r) n 



and N n -i = N n ~ t n - rj n l r 



In 

For k < n — 1, if Nk < tk we put eu = r\k = and A^-i = Nk- Otherwise we 
proceed as in the initial construction, i.e. put e/. = 1, rjk — [ Nk [~ tk J and A^_i = 
Nk — tk — rjklk- If Nk—i — we terminate and put ej = r\% = for all i < k. Since 
tk+i = tk + (L + 5 — l)lk and lk+i >tk>h for all k > 0, this algorithm yields a 
representation of N in the form 

n 

(5) N = Y J (eiU+r)ih), 

i=0 

where ej € {0,1}, < n < L + S — 2, and e; = implies 77^ = 0. Furthermore, since 
tk + (L — l)lk = lk+i, it follows that rji > L — 1 implies = 0. 
Note that the representation (O is not unique. Consider e.g. the case L = 2, 5 = 1; 
then we have t2 + h = 12 = ta + ii + Zi. However, for the rest of this paper speaking 
of a representation we will always mean the representation whose coefficients and 
r\i were constructed as explained in the algorithm above (i.e. in the above example 
the representation we choose is 12 = <2 + '2)- 

In order to establish unique "digit expansions" we use the following Lemma. 



s 



C. AISTLEITNER, M. HOFER, AND V. ZIEGLER 



Lemma 3. There is a bijection between positive integers and finite sequences of 
the form 

V = ((e n ,Vn), • • ■ , (eo,?7o)) 

such that e, G {0, 1}, e n = 1, < f] < L+S— 1, £j = implies rji — and rji > L — 1 
implies e i+1 = 0. TTiis bijection is given by 

a 
i=0 

and its inverse 

*- 1 (7V)-((e„,7 ] „),...,(e ,ry )), 
where the £j and 77, are computed by the algorithm described above. 

Proof. Let TV > be an integer. Note that \E' — 1 is injective since by construction 
we have Nk < ifc+i and therefore in every step of the algorithm the pair (e^, T)k) is 
uniquely determined. 

It remains to prove that 1 i r_1 (\f'('Z?)) = T). We prove this by induction on the 
length n+1 of V. The case n = is trivial, since ^(D) = £0 + ^0 < ti, and applying 
yields indeed T>. Let us assume that the algorithm yields the correct sequence 
^>(T>) for all sequences T> of length < n. In particular, this implies ^(D) < t n for 
all T> of length < n. Assume now that V is of length n+1. 

First, we prove that N = ^>(T>) < t n +i. Since e n — 1 we have T] n -i < L — 2 and 
by induction we know that for a sequence £ of length n starting with (1,L — 2) 
wehave *(£) < t n -i+{L-l)l n -i = l n . Hence * (V) < t n +(L+S-2)l n +l n = t n+1 . 

Now let £ be the sequence of length n induced by T> by deleting the entry 
(tmVn)- Note that e„_i might be zero and thus £ is not a valid output of the 
above algorithm. If ej = for all i < n — 1, then the proof is trivial. Assume now 
that at least one > for i < n — 1 . Then we can write £ = (£\ ,£2), where 
£\ = ((e n -ij Vn-i),..., (e n _ fe ,7? n _fe)) and (e^) = (0, 0) for n - fc < i < n - 1 
and £ 2 = ((e„-fc-i, r/„_ fc _i), . . . , (e ,??o)) with e„_fe_i = 1. We obtain that N = 
^(V) = t n + r] n l n + ^(£2) and since ^(£2) < In, we know that r\ n is the same 
integer which we obtain by applying our algorithm to N. In particular, we have 

= ((l,^),^,*- 1 ^))) = ((l, r]n ),£)=V. 

□ 

Using this digit expansions we are able to prove arithmetic properties of LS- 
sequences. We start with a lemma which provides conditions under which a point 
£l s nes m a cer tain elementary interval. 

Lemma 4. Let N be an integer with representation given in ([5]). Then 

n 

(6) = min { X ' e * + + m ax{e, + — L, 0}) . 

i=0 
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Moreover let 1% = [££ s>££ s + fi k )> with x < Ik, be an elementary interval. 
Then ^ s £ 1^ if and only if 

fc-i 

x = /Xejtj + f]tk) 

i=0 

is the truncated representation of N. 

In addition let Ai k \N) — (j{m : m < N,£™ s £ ii^} and assume that 

n 

N = x + sXejtj + rjik). 

i—k 

Then 

n—k 

Ci+kti + 1] i+ kk) + 1. 

z=0 

Proof. We start with the proof of (j6|). Due to Definition |4] we have 

g iS = (3 n min{L, e„ + + /3 ,l+1 max{e„ + Vn - L, 0} + ^ s 

with 

n-i 

^ = 5^(64*4 + ^). 

Repeating this argument inductively we will end up in (J6)). 
Now let N be of the form ((SJ) and let 

k — 1 n — fe 

M = y^(e^» + and 7V 2 = ^(e i+fc ti + r) i+k k). 

i=0 i=0 

Then by (|6]) we have ^ s = + /3 k £ L 2 s . Since the LS-sequences only take values 
in the interval [0,1) and since two distinct points of an LS-sequence with index < Ik 
differ at least by j3 k this implies the second statement of the lemma. Moreover, by 
assumption N% = x < Ik and therefore N2 can take all integer values, since we have 
no restriction on the digits tk and rjk (see Lemma [3]). Thus, -/V2 counts all points 

Q s in the interval 1^ with x < i < N . Since £ £ s £ Ix is the first point that 

(k) 

hits the interval I x , the last statement of the lemma is established. □ 

Note that LemmaS]would not hold in case of < x < tk, since this would imply 
£fe = (see Lemma and we would have serious restrictions for the digits of N 2 . 
Thus N2 could not take all integer values. 

Since a crucial point of the proof of Theorem Q] is to have a precise knowledge of 
(N) the next lemma gives a further method to describe this quantity. 

Lemma 5. Assume that x < Ik- If £l s ^ > then there exist integers A, B such 
that N = x + At k + Bl k and A {k) (N) — 1 + A + B. 
Proof. By Lemma 2] we know that N is of the form 

n 

N = x + ^2(e i t i + rnl i ). 

i—k 
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Now by the recurrence relations for tk and Ik and Lemma [T] we deduce that there 
are integers A, B such that N = x + Atk + Blk- 

Since s € 1% the proof of the lemma is complete, if we can prove that there 
exist integers A, B such that 

n n — k 

y^X^i + Vih) = At k + Bl k and ^ (ti+kU + r/ l+k h) = A + B, 

i—k i—0 

where e^, r/i and n > k are arbitrary non- negative integers. We prove this assertion 
by induction on n — k. The case n = k is trivial. We postpone checking the case 
n = k + 1 to the end of the proof. Now, let us assume that n — k > 2. Using the 
recurrence relations ([1]) and ([2]) for tk and Ik respectively we have 

n n—1 



y^X^j + Vih) = + Vih) 



and 



with 



i—k 



n—k n—k—1 



^(ei+fc^i + Vi+kk) — ^ (?t+fc*» + Vi+kk), 



ii = ej, rji — r\i for i = k, fc + 1, . . . , n — 3 and 

£n-2 = £n-2 + "Sen, ?n-l = Cri-1 + Le n , 
fjn-2 = r/n-2 + Sl] n , f) n -i = 7/ n _i + Ll] n . 

The new representations have fewer summands and by induction hypotheses we 
find appropriate integers A and B. 

It remains to check the case n — k = 1 . Using Lemma [1] we get 

^ktk+Vkh + £fc+l£fc+l + Vk+lh+1 

=tktk + Vkh + efc+iOfc + h(L + S - 1)) + %+i(t fe + l k {L - 1)) 



=t k (e k + fife+i + r] k+ i) +h (Vk + £fc+i(*i - 1) + Vk+i{h ~ 1)) ■ 

But, now 

A + B = 6 k + Efe+l + T] k+1 +T] k + e k+1 (t 1 - 1) + 7] k +\{h - 1)) 

= £fc<o + ??fe^o + (k+ih + Vk+ih 
and we have found appropriate integers A and B. □ 

The next lemma is related to the discrepancy of one-dimensional LS-sequences. 
In particular we are interested in an accurate formula for -^jj— -, where £^ s £ 1^ . 

Lemma 6. Assume that N has a representation of the form ([5]), and assume that 
s G /i*' 1 . TTien we /ia?je 

^ fc) (AQ _ flfe , gfl - (-sp) k ) + 1 - xp k 

N ~ P + TV 
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where 

n 

R = J2^iri + m^i)(SPy- k , 

i—k 

which can be estimated by 

\R\ < max{\n\, \ n + (L + S- 2)A 1 |} L^L_ . 

if S/3 ^ 1 and 

\R\ < max{|n|, |n + (L + S - 2)Ai|} max{n - k + 1, 0} 

if S/3 = 1. 

Proof. Using our assumptions and Lemma 2] we can calculate the exact values of 
A (k) (N) and N. In fact, we have 

n n 

N = x + 53(eitj + rjik) and A^(N) = Efe**-* + ViU-k) + 1- 

i— z— 

This yields 

^b fc) (AQ JL^ki^i-k + mii-k) + 1 

N J^^ki^i + Vih) + x 



i—k 



i—k 



k , ^R+l , Er= fc (^r + ^A )r' i+;c -^ fe 



N N 
(p- k Etki^o + Vi*o)p- i+k + + x) p k 

N 

k R{l-{-Sp) k ) + l-xf3 k 
1 N 
Thus it remains to estimate R. Note that t\ < and 

M + 77iAi| < max{|n|, |n + (L + S - 2)Ai|} 

for i = fc, . . . , n. Hence to complete the proof of Lemma[6]we only have to compute 
the geometric sum 

n— k 

El-* 

i=0 

We have to take absolute values in order to estimate R. □ 

Note that Lemma [6] gives a constant bound for \R\ for n — >■ oo if and only if 
Sj3 < 1, which is equivalent to L > S 1 — 1. This is exactly the case when we have a 
one-dimensional low-discrepancy LS-sequence. 
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3. Proofs of main results 

We start with the proof of Theorem [T] According to Definition [6l we define the 
recurrences 

f(i) - T ./W i O.+W +W - 1 +(*) _ r. 4. c. 

/(0 _ t ./W i o./(0 /(*) - i /(*) - r . 

<,(«) - r .,(*) i e. s W ,(*) - n s (l) - <?■ 

s n — ^«* n _i T "->j* n _2; *0 — u ' A l — °») 

which correspond to the number of intervals, long intervals and short intervals after 
n refinement steps in the i-th component of the two-dimensional LS-scquence, for 
i = 1,2. 

Now let k and m be integers satisfying the assumptions of Theorem Q] and as- 
sume that the integers k and rh are "large" (a precise condition will be given later). 
Furthermore choose two integers x\ < Zj^ and x 2 < Zm such that X\ ^ x 2 . In 

the sequel we will consider the intervals I = 1^ x 1^ and I = Ix^ +k ^ x I^ +m \ 
We want to prove that no point of the two-dimensional LS-sequence is contained 
in the interval I. Note that I C /, and consequently a point can only be contained 
in I if it is also contained in /. Now let N be given, and assume that £g <E I. We 
will also assume that £ I, and show that this leads to a contraction (provided 
k and to are sufficiently large) Consequently, no point of the sequence (£g)neN can 
be contained in /. 

Due to LemmaO s £ 1^ implies that there exist integers A±, B\ such that 
^ = 2^1+ A\v^ + Bil^\ Similarly, 5a £ iffi implies the existence of integers 
A 2l B 2 such that N = x 2 + A 2 t$ + B 2 l$ . Thus 

(7) x x + Axt^ + BilP =x 2 + A 2 t$ + B 2 lg>. 

Moreover, choosing A\ and B\ according to Lemma [5] we know that there are 
exactly 1 + A\ -\ 
Lemma [5] yields 



(k) 

exactly 1 + A\ + B\ points with index < N lying in the interval Ix x ■ Therefore 



Ax+Bx + 1 _ Rk R x {\ - (-StPi)") + 1 - XxPt 



xi + Axtr + Bit 



k T ^^k 

Multiplying both sides with N = x\ + Ait£' + Bilj}' and solving for B\ we obtain 

i - i - 

A similar argument for the second component yields 

(9) B 2 = A 2 — 1 — . 

1 pm 1 j\^J om 

1 — l m p 2 1 i k p 2 
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Now we resubstitute equations ((8]) and dH) into (0) and obtain 

1.(1) D ;(l)n _ f_Q. «2\fe 



1 - 1 - 



x 2 + A 2 



1 — R m 1 — ;( 2 ) «m 



Now let us investigate the quantities prj^fs an d i-h/f k 
Lemma 7. FFe /ictue 



and 

lk{ \ { ff h 2)k) =^- k + hs^. 

1 - l k p h Ai 

Proof. Using the explicit formulas ([1]) and ([2]) for the recurrences tt and Ik, respec- 
tively, we obtain 

h - h Jtq - X Q )f3- k + (n - A!)(-S/3) fc = Tjfeg^ - (-^ fc ) 
l-^ fc l-p k (X p- k + Xi(-Sp) k ) A!(l-(-5/? 2 ) fc ) 

VL^+4g o_fc 1 - (-S/3 2 ) k = o-fc-1 

Ai ' P ' 1 - (-5/3 2 ) fc P 
which proves the first part of Lemma [7] 



Note that 



-L+VL 2 +4,S 



Ai _ -L + VL 2 + 45 = 



25 



V£ 2 +4S V£ 2 +4S 
Let us now prove the second statement of the lemma. Again we use the explicit 
formulas (JTJ) and ([2]) and obtain 

l k (l - (S(3 2 ) k ) l k (l - (-5/3 2 ) fc ) 

1 - Zfcy5 fc l-/3 fc (A /?- fe + A 1 (-5/3) fe ) 

I fc (l - (S(3 2 ) k ) 



Mi-(-W) 

+ (-5 / 9) fc 



□ 

Continuing the proof of Theorem[Tl let us insert the explicit formulas of Lemma[7] 
into mOJl. We get 

"PI 



(11) A 1 =A 2 g- I + (x 2 ~x 1 )p k+1 +R, 



qm-{-l 

where 



v( 2 ) \ flfc+1 /\(!) 



(12) R = R 2 p 2 ( ^ + (-$ 2 $) m J A_ ft ^ + (-5i^) fc ) 
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ok + l 

Now by assumption we have = „ f° r some coprime positive integers p, q 

and (fTTj) can be written as 

(13) 2(qA 1 - P A 2 ) = 2{x 2 - x 1 )qf3 k+1 + 2qR, 

Obviously the left side of (fT3|) is an even integer. In order to prove Theorem Q] we 
want to show that the right side is not an even integer if £g G I and k and fa are 
sufficiently large. By our assumptions on x\,x-i and the parameters L\,S\,Li and 
S 2 we have 

(14) ||2 g (.T 2 - a;i )/3 1 fe+1 || odd = e< 1, 

where || • \\ dd denotes the distance to the nearest odd integer. Note that /3^ +1 is 
irrational due to Lemma [5] and therefore indeed e < 1 . 

Now Lemma [5] tells us that the assumption G I yields 

n 

i—k-\-k 

and we have 

1*1 < max{|r«|, |r« + {la + * - 2)A«|}t^£. 

A similar inequality also holds for |i? 2 |. Now we choose k sufficiently large such 
that 

\Si(3if < m if-^ 7fl-(l - S1P1) 

Idl^maxllr^l, + {L x + S x - 2)A«|} 

in case of |ci| ^ and k = otherwise. Similarly we choose m sufficiently large 
such that 

< ^ ^ 7=T-(1 - S 2 fc) 

|c 2 |4gmax{|r 1 (2) |, |r x (2) + (L 2 + S 2 - 2)X^\} 
in case of |c 2 | 7^ and rh = otherwise. With this choice of k and fa we obtain 

I^K 1 ' 



2q ' 

which together with (fT4|) implies that the right side of (|13|) is an odd integer plus 
something < 1 in absolute values. Consequently (fTTj) does not have any solution 
Ai,A 2 . However, since (fTT|) followed from our assumptions, this means that we 
have a contradiction. Consequently it is not possible that £g G / for any N, which 
means that I does not contain any point of (£g)neN- This completes the proof of 
Theorem [TJ 

We continue with the proof of Theorem [5] Let b — gcd(L 1; Si,L 2 , S 2 ) > 1 and 
let t [ n ] , l ( n 1} , t ( n ] , l£ ] be defined as in the proof of Theorem [TJ Note that b divides 
t { n\ln\tn\ln ] for all n > 1. Now consider the interval I = [0, x [f3 2 ,2f3 2 ). 
Note that 2j3 2 < 1 since by gcd(Li, Si, L 2 , S 2 ) > 1 it follows that L 2 > 2. It 
follows from Lemma [5] that we can write every N± for which £z/s € [0, /3i) as 
^1=0 + A^P + Bil^ and every N 2 for which ^.s, G [/?2,2/3 2 ) as N 2 = 

(2) (2) 

1 + A 2 t\ + for appropriate integers Ai, Bi, A 2l B 2 . Hence we obtain for 
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every N±,N2 that N% = mod b and N2 = 1 mod 6, respectively. Consequently 
no point of the two-dimensional LS-sequence can be contained in /. This proves 
Theorem [51 
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